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\ Abstract 

' We analyze some consequences of grand unification of the third-generation 

Yukawa couplings, in the context of the minimal supersymmetric standard 
model. We address two issues: the prediction of the top quark mass, and 
the generation of the top-bottom mass hierarchy through a hierarchy of Higgs 
vacuum expectation values. The top mass is strongly dependent on a certain 
ratio of superpartner masses. And the VEV hierarchy always entails some 
tuning of the GUT-scale parameters. We study the RG equations and their 
semi-analytic solutions, which exhibit several interesting features, such as a fo- 
cusing effect for a large Yukawa coupling in the limit of certain symmetries and 
a correlation between the A terms (which contribute to b — > 57) and the gaug- 
ino masses. This study shows that non-universal soft-SUSY-breaking masses 
are favored (in particular for splitting the Higgs doublets via D-terms and for 
allowing more natural scenarios of symmetry breaking), and hints at features 
desired in Yukawa-unified models. Several phenomenological implications are 
also revealed. 
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1 Introduction 



There is strong evidence to suggest that the three gauge couplings of the strong, elec- 
tromagnetic and weak interactions are unified at a high energy scale in a single gauge 
interaction based on a simple group, such as SU(5) or SO(10), as long as the desert 
below the unification scale is described by the minimal supersymmetric extension of 
the standard model (the MSSM). Furthermore, the combination of supersymmetry 
(SUSY) and grand unification yields models with numerous attractive features: the 
embedding of the standard-model matter multiplets into a few irreducible representa- 
tions of the GUT group, the technically natural preservation of a hierarchy between 
the weak and GUT scales, a longer proton lifetime to allow agreement with current 
experimental lower bounds, the correct prediction of the ratio of b quark to r lepton 
masses, simple ansatze for the remaining fermion masses, and a picturesque scenario 
for radiative electroweak symmetry breaking. We have chosen, therefore, to look 
beyond the gauge unification prediction of the weak mixing angle and examine the 
unification of third-family Yukawa couplings Jl] , for the most part within its natural 
context of SO(10) unification. By Yukawa couplings we mean the couplings of the 
top, bottom and tau to the Higgs doublets which generate their masses when elec- 
troweak symmetry is broken. The third family is singled out because of its relatively 
large Yukawa couplings: it seems reasonable to suppose that they arise at tree-level 
from the simplest interactions, while the masses and mixings of the other generations 
require more complex, perhaps also higher-order and certainly very model-dependent 
structures. Our study also applies more generally to scenarios in which these Yukawa 
couplings are unified but not in the context of an SO (10) GUT, or even models in 
which the Yukawas are only comparable near the GUT or Planck scales. We begin 
by considering the most immediate prediction of this Yukawa unification, namely the 
top quark mass ||. However, we are quickly led to consider in some depth the more 
general question of how the top-bottom mass hierarchy could be generated in the 
MSSM, and how this hierarchy depends on the initial conditions of the renormaliza- 
tion group (RG) evolution at the GUT scale M. We will conclude with a discussion 
of how natural (or unnatural!) such a hierarchy seems in this context, what its other 
phenomenological predictions might be, and how one could hope to improve the the- 
oretical picture or obtain experimental corroboration. 

One consequence of Yukawa unification is immediate, and independent of other 
assumptions except for the qualitative nature of the RG evolution equations of the 
MSSM. Since the Yukawa couplings of the top and bottom quarks (and the r lepton) 
are always comparable, the large ratio of the top mass versus the bottom (or tau) 
mass must be due to a large ratio of the Higgs vacuum expectation values (VEVs) 
which give rise to their masses. Namely, since the up-type and down-type matter 
fermion masses arise from couplings to the up-type and down-type Higgs multiplets 
(Hu t £>), respectively, the large ratio m t /m hT = (\t v u)/{\,T v D) is n °t a consequence 
of a large ratio of Yukawas Xt/Xb,r but rather of large vjj/vd = tan/3. Thus Yukawa 
unification generically implies tan/3 ~ O(50). Further assumptions are necessary to 
make any precise predictions. We will assume the following three throughout most 
of this work, although we will point out those conclusions which are more general: 
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(I) The masses of the third generation, m t , mb and m T , originate from renor- 
malizable Yukawa couplings of the form I63 O 16 3 in a supersymmetric GUT 
with a gauge group containing (the conventional) SO(10); I63 denotes the 16- 
dimensional spinor representation of SO(10) containing the third-generation 
standard-model fermions (plus the right-handed neutrino which we assume to 
be superheavy) and their superpartners. 

(II) The evolution of the gauge and Yukawa couplings in the effective theory beneath 
the SO(10) breaking scale is described by the RG equations of the MSSM. 

(Ill) The two Higgs doublets lie predominantly in a single irreducible multiplet of 
SO(10). 

The first and third assumptions serve to define what we mean by Yukawa unifica- 
tion, while the second allows us to relate this GUT-scale unification to weak-scale 
observables. From (I) and (III) it follows that the third-generation Yukawas must 
arise from either a 16 3 10 g 16 ? or a 16 ? 126 g 16 ? interaction with an SO(10) Higgs 
multiplet. The latter leads to the boundary conditions 3A^ = 3A^ = \^ = \ G at 
the GUT scale, but the resulting ratio of m h /m T at low energies is far too low to 
be consistent with experiment (at least within the perturbative regime, and unless 
very large threshold corrections to the h mass arise at low energies 0]). Thus we are 
restricted to using the 10 H , and hence the boundary condition 



With this boundary condition, and using the unification of gauge couplings to fix 
the unification scale and the gauge coupling at that scale, we can now evolve the 
Yukawa couplings down to the weak scale for any given value of Xg- The idea is that 
the three observable masses m t , m>b and m T are functions of the four parameters A t , 
Xb, A T and tan /3, and A^& jT are in turn determined by the unification in terms of \q 
and the GUT scale Mq- Since the latter is already known from gauge unification, 
we are left with three observable masses as functions of only two parameters, Xg and 
tan/?. Thus we use two observables, nib and m T to fix Xg and tan/5, and thereby 
predict the third observable m t . A detailed analysis of the RG evolution, and the 
consequent predictions, has already been presented 0J. The results, namely the values 
of X t ,b,r and the ratio R = rrib/m T = Xb/X T all at the weak scale, are plotted in Fig. 1 
as functions of Xg- (These curves actually use 2-loop RG evolution, but at this 
point the difference between 1- and 2-loop equations is not important. For the final 
predictions of m t we use the full 2-loop evolution and 1-loop matching conditions.) 
Evidently, larger Xg values correspond to a heavy top and to a smaller R ratio. The 
experimental value -R exp t is found 0J from the QCD sum rules value and is evolved to 
the weak scale using 2-loop QCD running. We find, allowing for a s to vary between 
roughly 0.11 and 0.12, the shaded range shown in Fig. 1. Thus, in the absence of any 
large corrections in the matching between the R evolved down from the GUT scale 
and the -R exp t in the standard model, we find A^ > 0.75, which implies a heavy top. 
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2 The top and bottom masses 



For a precise prediction, 2-loop RG equations must be used along with 1-loop match- 
ing functions. These matching functions include logarithmic corrections from infinite 
counterterms as well as nonlogarithmic contributions from finite graphs. The former 
are given elsewhere 0; they are generally quite small, and invariably increase the top 
mass as the superpartner masses increase. The latter are more interesting, since they 
can be very large || The dominant corrections arise from the graphs of Fig. 2, 
which match the value of the b mass as evolved down from the GUT scale to the value 
in the low-energy theory. Typically the gluino graph dominates, yielding a corrected 
value m b = X^Vu + 8m b where vjj = 174 GeV, 

5m b 8 2 tan (3 rrigfj, 

= o 93 77 9 T~ ' v ) 

m b 3 I67H m^j 

and trig is the gluino mass while m e g- is the mass of the heaviest superpartner in the 
loop (more exact expressions may be found in our previous work). The important 
observation here is that the bottom mass gets a large contribution from the up- 
type Higgs at 1-loop order, whereas its tree-level mass was small due to the small 
VEV of the down-type Higgs. In the usual scenario with small tan/3, the bottom 
was light because its Yukawa coupling was small, or in other words because it was 
protected by an approximate chiral symmetry. Thus any higher-order corrections 
would also be suppressed by this approximate symmetry. In the large tan j3 scenario 
these corrections are not suppressed: there is an enhancement of Vu/vn = tan/?, 
which overcomes the usual gf/167r 2 loop factor to give a correction of order 1 to the 
b mass — at least if rrigfi ~ m ls- Phenomeno logically, the result is that m b cannot 
be predicted with any certainty unless we know something about the superspectrum, 
and if m b is uncertain then so is the top mass prediction. 

In Fig. 3 we show the top pole mass prediction, now to full 2-loop order, as a 
function of the MS running parameter m b (m b ) in two cases. The top curves and the 
higher horizontal axis correspond to a hierarchical spectrum, in which the squarks are 
heavy whereas the \x parameter and gaugino masses are light. Then the corrections 
5m b are small and the top mass is predicted to be above 180 GeV or so. The bottom 
curves and the lower horizontal axis correspond to a roughly degenerate spectrum in 
which the corrections to the b mass are (9(25%) and negative (i.e. R should be lowered 
by 25% in Fig. 1 before matching to the experimental value). Now the top can be 
significantly lighter. In fact, this last argument can be turned around: if the threshold 
corrections are too large in magnitude and negative then the top mass prediction will 
be below the experimental lower bound, while if they are large and positive then no 
value of Xg will allow agreement with -R ex pt- We thus find the following bounds on 
the superspectrum: if Xg is allowed to vary, then 

_0.37<^<0.08, (3) 
whereas for fixed Xg the appropriate limits can be read from Fig. 1. 
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3 Radiative bottom decay (I) 



Do we have any experimental information about these corrections? Of course we have 
no direct evidence of any of the superpartners, but we can appeal to their indirect 
appearance in loop diagrams. Consider again the diagrams of Fig. 2, but with one 
of the b quarks replaced by a strange quark using a flavor-changing vertex, and with 
a photon attached in all possible ways. We see that the same processes will lead 
to a contribution to the rare decay b — > 37, and with a similar enhancement of 
C(tan/3) over the usual MSSM scenario |2|, ||. The dependence of these diagrams 
on the superpartner masses is somewhat similar to that of Srrib, except that (a) 
now the operator is of higher dimension and so is suppressed by the mass of the 
heaviest superpartner, and (b) typically the higgsino-mediated diagram dominates. 
If the parameters appearing in this Higgsino diagram, namely /i, A (the trilinear soft 
SUSY-breaking parameter) and the squark masses, are all comparable and of order 
the Z mass, then this diagram gives a contribution to the amplitude for b — *> s'j many 
times bigger than the standard model or the usual MSSM amplitudes, and is clearly 
ruled out by the CLEO limit ||. To restore agreement with experiment, either the 
overall superpartner mass scale must be raised far above the electroweak scale, or else 
fi or A or both must be suppressed relative to the squark mass in the loop. More 
quantitatively, we find that either the masses must all be raised to at least C(TeV), or 
else if both [i and A are near the Z mass then the squarks must be above either ~ 400 
GeV or ~ 700 GeV, depending on whether these diagrams interfere destructively or 
constructively with the 2-Higgs standard model amplitudes. In any case, however, 
these restrictions do not yet tell us anything about the combination mgii/m 2 s which 
appears in the link between these two will be forged below, when we study the 
evolution of the entire set of MSSM parameters. 



4 Electroweak symmetry breaking 

We have assumed in the above that the top-bottom mass hierarchy would arise from 
a hierarchy of VEVs in the Higgs spectrum, namely Vujvo = tan/3 ~ (9(50). We 
now consider how such a hierarchy could be generated in the MSSM. Clearly this 
question is of interest for any model in which the Yukawas themselves do not supply 
a sufficient hierarchy to explain the large ratio of third-generation quark masses, 
not just for the equal- Yukawas case to which we have specialized. In studying this 
question, however, we will pay attention to which spectra are favored by large tan f3 
scenarios, and whether with such spectra and our unification assumptions we can pin 
down the top mass prediction. 

The scalar potential of the neutral Higgs bosons which leads to electroweak sym- 
metry breaking is given by 

V Q = mljlHul 2 + m 2 D \H D \ 2 + ^B(HuH D + h.c.) + (quartic terms) (4) 

where mfj D = fifj D + fi 2 contain the soft-breaking masses and the /1 parameter in the 
superpotential, B is a soft-breaking mass parameter, and the quartic terms arise from 
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D-terms and so are given by gauge couplings. We will evolve these parameters from 
the GUT to the weak scale using the 1-loop RG equations of the MSSM, stopping 
the evolution at some typical scale (of order the squark masses) which minimizes 
the effects of higher-order corrections. The subscript indicates that we will restrict 
our attention to this (RG-improved) tree-level potential rather than calculate the full 
1-loop effective potential or, better yet, explicitly integrate out massive particles and 
consider full 1-loop matching conditions. We expect that our qualitative discussion 
of the radiative symmetry breaking will not be jeopardized by this simplification. 
That is, a more complete calculation will change the numerical values of the GUT 
parameters needed for correctly breaking the symmetry, but will not significantly 
alter the size of the domains in parameter space where such breaking is achieved. 
The conditions for this breaking are well-known: 

ml + ml >2\/iB\ (5) 

ensures that the potential is bounded from below, and 

m 2 um 2 D < fi 2 B 2 (6) 

guarantees the existence of a minimum away from the origin and so breaks the sym- 
metry. In practice, since \fiB\ will always be much less than or at most comparable 
to \m\j\ and |m|,|, we can reduce these requirements to m 2 A = mfj + m 2 D > (using 
the expression for the pseudoscalar Higgs mass) and mfj < (noting that large tan j3 
means that the up-type Higgs gets the large VEV). 

In the usual — and very attractive — scenario of radiative breaking, the two mass 
parameters start out at the GUT scale with a universal positive value: mfj = m 2 D = 
Mfj + /i 2 , where M H is the soft-breaking mass of the 10# of Higgs in SO (10), or of the 
5 H and 5^ in SU(5). Thus the symmetry is not broken at that scale. However, in the 
RG evolution to the electroweak scale, the large Yukawa coupling of the top quark to 
Hu, which gives the top its mass, also drives the mass-squared parameter fifj of Hjj 
negative (with the help of the QCD coupling), while the absence of a large Yukawa in 
the down sector keeps the mass-squared of H D positive. In fact, conditions (j5|) and 
(0) are easily satisfied for a large range of initial conditions if \q ~ 0(1), resulting 
in a very natural picture of radiative symmetry breaking. This picture is essentially 
lost in the large tan f3 scenario, for the following two reasons: 

1. Since both Yukawas are comparable [and in fact initially equal in the SO(10) 
case], the two Higgs doublets tend to run in the same way, so either both stay 
positive at the electroweak scale and the symmetry does not break, or both 
become negative and the potential becomes unbounded from below (a situation 
which breaks the symmetry but in a Coleman- Weinberg fashion, yielding an 
"electroweak scale" orders of magnitude higher than the SUSY-breaking scale). 
The effects which differentiate the evolution of the two Higgs doublets, namely 
hypercharge and the absence of a right-handed neutrino, are small and a poor 
replacement for the usual At ^> A& splitting. Interestingly, an 0(1) splitting 
between A f and \ T is still of little use since it is quickly diminished by the 
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fixed-point behavior of these couplings. (Some of these observation have been 
previously made by T. Banks 0.) 

2. Even when electroweak symmetry is broken, a large hierarchy of VEVs must be 
generated between the two similarly-evolving Higgs doublets. Minimizing the 
potential Vq when tan j3 > 1 yields 



and 



- 2ml = m% (7) 

tan (3 m\ + m 2 D m\ 

The first equation sets the scale, but from the second equation we see that a 
large hierarchy in VEVs requires the large hierarchy fiB <C fnfj + m 2 D . This, as 
we show below, implies a degree of fine-tuning between some parameters in the 
Lagrangian. 



5 Solutions of the RG equations (I) 

Before analyzing the implications of these two criticisms, we present the 1-loop so- 
lutions || of the RG equations for the MSSM mass parameters, integrated between 
M G = 3 x 10 16 GeV and a typical squark mass of 300 GeV. (None of our results is 
sensitive to the exact values of these starting and stopping scales.) The solutions 
depend on the dimensionless initial values of the gauge and Yukawa couplings «g and 
X G and on the dimensionful GUT-scale parameters M sq , M H , fj, G , Mi/ 2 , A G , B G and 
M X - M sq and M H are the soft-breaking masses of the 10# and the I63 respectively, 
and Mx will be explained below; we note for now that it vanishes for universal soft- 
breaking masses. Whenever possible, capital letters denote values at the GUT scale. 
The RG equations themselves are well-known and will not be presented here. These 
equations dictate that the low-energy values of the various mass parameters depend 
very simply on the dimensionful initial values, with coefficients that depend only on 
the dimensionless ones. Since a G is known from gauge unification, these coefficients 
depend only on \ G . For the representative value \ G = 1, the solutions are: 

2ml = - 5 - 1 M \I2 + 1 - 2 Mf -1.6Mg 2 + 2/J 2 -3.8 Mf (9) 

m\ = -4.9Mg 2 2 + l.lMf 2 -1.7M s c f + 2/i 2 + .01M<f (10) 

m 2 Q = +4.6Mg 2 2 -.25M<f + .51MJ 2 +1.0Mf 2 (11) 

m\ = +4. 1 Mf /2 2 -. 27 M<f + . 46 Mg 2 + .85 Mf 2 (12) 

m 2 h = +4.2 Mg 2 2 - .23 Mf + .55 Mf q 2 -2.9M| 2 (13) 

m 2 L = +.53Mg 2 2 -.12Mf 2 + .77M^ 2 -3.1Mjf (14) 
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2 - +.15Mg 2 -.23M| +. 55 M a <f +l.2Mf (15) 
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G 

1/ 



A t = +.09A G + 1.8Mg 2 (16) 
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A b = +.07A G + 1.9Mf /2 (17) 

A T = +.20A G -.17Mg 2 (18) 

B = -.86 A G - l.lMg 2 + 1.05 G (19) 

/i = -44 /i G (20) 

Here Q and L are the squark and slepton doublets, respectively, and t, b and r are 
the SU(2)-singlet squarks and sleptons. For clarity of presentation, we have dropped 
from the first 7 expressions above the terms proportional to A 2 G and to AgMi/ 2i 
since their coefficients are small enough [0(0.01 — 0.1)] and exhibit sufficiently small 
custodial-SU(2) breaking to be negligible for purposes of symmetry-breaking, at least 
if Aq is not very much larger than the other mass parameters. We have also used 
the low-energy value of /i in Eqs. (|^- [1Q|) . These solutions are useful references for the 
discussions below. These solutions can also be combined [Q] with the more analytic 
approach briefly described in Eqs. (|25|) to give a complete, semi-analytic solution to 
the RG equations when custodial SU(2) is an approximate symmetry. 



6 Obtaining a hierarchy of VEVs 

To understand the second criticism above, let us examine in more detail how Eq. (^) 
may be satisfied. We concentrate for the moment on six relevant electroweak-scale 
parameters: the up- and down-type Higgs masses and m 2 Dl a typical squark mass 
Mq, the B parameter in the scalar potential, a gaugino mass (specifically the wino 
mass) mi/2, an d the /i parameter. This last one may be set to zero by imposing a 
Peccei-Quinn symmetry on the Lagrangian, so the size of /i measures the breaking of 
this VQ symmetry; consequently, /i is multiplicatively renormalized. The previous 
two, B and m^, along with the A parameter, transform in the same way under a 
continuous 1Z symmetry (so they enter into each other's RG equations) and may be 
made arbitrarily small by imposing this 1Z symmetry. With these two symmetries in 
mind ||, we consider three possible spectra having splittings which lead to a large 
tan/? according to Eq. (||): 



mass: 


scenario A: 


scenario B: 


scenario C: 


7m z 








m z 


m;y m D 777o 


muin D mQ m 1 / 2 fJ> 


Vflu BlTLi/2^ 




Bm 1/2 /J, 










B 





Of course there are many other ways to split the parameters and obtain the correct 
hierarchy, but these will suffice to demonstrate the fine-tuning involved in the split- 
tings. The value of tan ft is determined directly only by mu, m D , \i and B; we include 
m and m 1 / 2 to illustrate the symmetries. Scenario A involves no tuning at all (at 
this stage of the analysis): the only hierarchies present are those enforced by the two 
symmetries. However, as also pointed out by Nelson and Randall 0, such a scenario 
is ruled out for large tan ft since it would imply a light chargino, in disagreement with 
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bounds from LEP ||. In fact, both \x and myj must be comparable to or above the Z 
mass to satisfy this bound. Therefore we must widely split some parameters without 
a symmetry justification, and this will entail fine-tuning. For example, in scenario B 
all parameters are kept at the Z mass but the B has been tuned to be light (namely, 
its initial value at the GUT scale is chosen to almost completely cancel the contribu- 
tions induced by A and mu2 through the RG evolution). Alternatively, in scenario C, 
mjj is adjusted to end up much below the other scalar masses (yielding m\ ~ 50m 2 z ) 
while the other parameters are kept at the Z mass using the approximate symmetries. 
In these two scenarios, and in fact generically whenever \x > mz and mu2 > mz, the 
initial conditions at the GUT scale must be adjusted to at least a relative accuracy of 
1/ tan/5 to obtain the necessary hierarchy of VEVs. We should point out, however, 
that such a tuning is no worse than the one which would be needed in the small tan f3 
case if the squarks were experimentally determined to be above 700 GeV or so. 



7 Splitting the Higgs doublets 

We return now to the first criticism above, and address the splitting between the two 
Higgs doublets. Recall that after running we need Imfj < while mfj + m 2 D > 0. 
However, the two masses evolve almost in parallel, since custodial symmetry breaking 
effects, namely hypercharge and the absence of vr, are small. Thus, if at the GUT 
scale the mass parameters are custodial-SU(2) symmetric, the splitting of the two 
Higgs masses at the weak scale is small relative to a typical SUSY mass Mg at the 
GUT scale: m 2 D — mfj = e c M| ("c" for custodial). Putting these together, we learn 
that only within a window of size ~ e c in the GUT-scale parameter space can we 
simultaneously satisfy mfj < and m 2 A > 0; if they are satisfied, then m 2 z = — 2m 2 J < 
e c M| and m 2 A = mfj + m 2 D < e c M|. In practice, this is usually accomplished |1[ 
using the gaugino mass as the largest mass parameter, so M s = M 1 / 2 > M sq>H : this 
is because, according to Eqs. (P|,|TU|), custodial breaking effects proportional to Mjy 2 
lower mfj with respect to m 2 Dl while those from the scalar masses M 2 qH act in the 
opposite way. Furthermore fi must also typically be 0{Mi/ 2 ) in order to keep m 2 A 
positive. Then, in addition to the 0(e c ) fine-tuning of the Z mass, the large values 
of mi/2 and /i mean that the B parameter must be adjusted beyond the 0(1/ tan/3) 
accuracy of the previous paragraph. To see this, we rewrite Eq. (||) in the form 

B _ 1 mfj + ml 
mx/2 tan/3 /itn^ 

which quantifies the needed suppression of the electroweak-scale value of B (achieved 
by fine-tuning its GUT-scale value) relative to the minimum value it would naturally 
have, namely the value ~ Mi/ 2 induced through the RG evolution. In the present 
case, using /z ~ Mi/ 2 ~ Ms we obtain B/mi/ 2 ~ (1/ tan/5) e c . 
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8 D-terms 



This highly unnatural state of affairs arises partly because of the degeneracy of the 
Higgs doublets and their subsequent parallel evolution. A possible remedy is actually 
generic in SO (10) unification, due to the rank of this group which exceeds by one 
the rank of SU(5) or the standard model. Thus we write SO(10) D SU(5) Cg> U(l)x, 
where U(l)x is proportional to 3(B — L) + 4T3R [the generator of baryon- minus 
lepton-number symmetry and a generator of SU(2)r] and couples to the scalar fields 
according to the following table: 



field: 


H v 


H D 


Q 


t 


b 


L 


T 




Ofiff) 


U(l)x charge: 


-2 


2 


1 


1 


-3 


-3 


1 


5 


-5 



The 16 H and 1Q H are examples of extra superheavy Higgs representations which are 
typically added in order to break this U(l)x (in this case by acquiring VEVs in 
the "z/r" direction) and reduce the rank of the group. As usual, the spontaneous 
breakdown of a U(l) leads to a VEV for its D-term, which can induce masses for 
the various fields which appear in this D-term. In particular, if we do not assume 
universal soft-breaking masses for all scalars, then the soft-breaking masses of the 16 H 
and 16 H need not be equal, and therefore their VEVs are also split, in proportion to 
their mass splitting. This splitting in turn generates a mass splitting in the low-energy 
MSSM Lagrangian through the cross-term: 

£ 3 \d\ = ~ ((\16 H \ 2 ) ~ (m H \ 2 ) + 2\Hu\ 2 ~ 2\H D \ 2 + . . .f . (22) 

In fact, this mechanism splits any fields which have different charges under U(l)x- 
Thus the boundary conditions for the scalar masses at the GUT scale become 

My = M 2 H + /1 2 -2M 2 X 

M 2 D = M 2 H + /1 2 + 2M 2 X (23) 

M 2 QAT = M 2 q + M\ 

M 2 L = M 2 q -3M 2 

where the capital letters on the left-hand side serve as reminders that these are the 
values at the GUT scale, and 

Ml = ^(M 2 6 - M^) (24) 

is a new soft-breaking mass parameter in the low-energy theory. With this mass 
we no longer need rely on large gaugino masses to split the Higgs doublets: they can 
start out being different, and thus even with parallel evolution the symmetry-breaking 
conditions ([§-§]) can apparently be satisfied. 

One problem with this mechanism is evident from the initial conditions in Eq. ([23]) : 
not just the Higgs doublets but also the squarks and sleptons are split, so an exces- 
sively large M\ could lower M 2 or M| sufficiently to make m 2 or m 2 L negative at 
the electroweak scale, thereby spontaneously breaking the strong or electromagnetic 
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gauge symmetries. If RG effects were irrelevant, namely for small Xq, then M sq 
could always be raised enough to prevent this without affecting mfj D . However, for 
\ G r^j (9(1) the squark masses strongly affect the evolution of the Higgs doublets 
[see Eqs. (ffllDDI, and only for very constrained ranges of the initial parameters can 
the electroweak symmetry, and only that symmetry, be spontaneously broken at a 
reasonable scale. In fact, as we show in brief below, there is a focusing effect that is 
inherent in the MSSM RG equations when At, ~ X t , and which inevitably requires an 
adjustment of the GUT-scale parameters beyond the 1/ tan/3 level derived above. We 
first show this behavior of the RG equations for completely general initial conditions, 
and then return to discuss the specific case of Eqs. fl23|) . 

9 Solutions of the RG equations (II) 

Consider the RG equations of the MSSM in the limit of exact VQ and TZ symmetries, 
in which /i = Mi/ 2 = A = B = at all scales. For future reference, we call this 
scenario the maximally symmetric case. This limit is interesting for two reasons: First, 
no large corrections arise to the b quark mass, and the R = m\,jm T prediction for all 
values of A^ ~ falls nicely within the range allowed by experiment (see Fig. 1); 
in other words, a heavy top quark near its fixed-point mass favors small bm^. Second, 
as we saw above, having a large fi and mi/2 calls for fine-tuning B (or some equivalent 
adjustment), so we'd like to explore the opposite limit to see whether a more natural 
scenario can be achieved. Of course, eventually we must relax this limit to agree with 
LEP bounds, but the qualitative behavior we shall discover will persist. If we further 
approximate A;, ~ A t = A and neglect the sleptonic contributions (thereby restoring 
custodial symmetry), the RG solutions simplify considerably. There are now five 
relevant parameters. In terms of their initial conditions at the GUT scale, M^, M|,, 
Mq, M| and M fc 2 , the solutions at the electroweak scale are: 




-f e x X - \e' x X' - I 
le x X + I 

^t\X — \l 

\t\X + \t'x x ' - 4I 
\t\X - \e' x X' - \I 



+ 




r 



(25) 



where 




(26) 



0.085 



(for Ag ^ 1) , 




(27) 



and 



X 



Ml + M 2 D + 2M% + Ml + Ml 
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X' 

I 
r 
i" 



Ml 

m 



Ml 
it + Ml) 
-Ml) 



u 



+ M? 



Ml 



- f(2Mj + M 2 t + Ml) 



(28) 



mi 



Ml) 



2M 2 - Ml - Ml 



(Note again the use of capital letters to denote GUT-scale initial parameters, and 
recall that U, D, Q, t and b refer to the up-type Higgs, the down-type Higgs, the 
SU(2)-doublet third-generation squarks, the SU(2)-singlet stop and the SU(2)-singlet 
sbottom, respectively.) 

Evidently, two linear combinations of masses, labeled by X and X' at the GUT 
scale, renormalize multiplicatively and exponentially contract at low energies for \q ~ 
(9(1). The three other linear combinations, /, /' and I", are invariant. The important 
observation here is that in the first contraction — the sum rule m\ + 2m 2 Q + m 2 + 
ml = e\X — the coefficient of every term is positive, while we already know that 
each mass-squared itself must be positive for a proper electroweak-breaking scenario. 
Therefore each term by itself must be small, less than t\X. This can only happen 
if the various combinations of invariants (and possibly also e' x X f ) in the expressions 
(p5|) for these terms are adjusted to be small relative to X. The exact constraints 
that follow from this requirement are given explicitly elsewhere ||. They are of the 
form J, I', I" < ma.x(e\X,e' x X'). We learn that, for Ac ~ C(i) where this focusing 
effect is important, any given model for the GUT-scale soft-breaking masses must 
either provide an explanation of why each invariant should be small relative to the 
sum X = Mfj + Ml + 2Mq + Ml + M 6 2 , or else that invariant must be tuned by hand 
to be small. We also learn that the conditions for successful symmetry-breaking are 
sensitive to any other small effects. One such effect is custodial symmetry violation, 
which is parametrized above by e c and results from hypercharge and A T (or the absence 
of z/r). In the running of the Yukawas, both of these cause X t to slightly exceed 
and thus drive mfj below m 2 D , as in the conventional scenario of electroweak symmetry 
breaking. In the running of the masses, the contributions of the r Yukawa have an 
opposite and numerically more relevant impact. Therefore the custodial-breaking 
effects make it harder to break the symmetry correctly — significantly harder in the 
specialized case discussed below. In any realistic scenario there are also contributions 
from the gauginos and /i, so in the end the sum rule, and therefore the general limit 
which must be set on J, I' and J", takes the form 



{l,l',l",m 2 A 



2m 2 Q + m\ + ml] < O 



max e\,e c 



Mr 



1/2 



Ml 



M 



(29) 



where Ms is the largest mass parameter in the initial conditions at the GUT scale. 

If we now return to the more specialized boundary conditions of Eq. (|23|), we find 
(after setting /i = 0): 



X 
X' 

I 



2M 2 H + m% 



l(2M 2 H -3Ml q ) 



(30) 
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V = -AM 2 X 
I" = +AMl 



For this choice of boundary conditions, keeping the invariants small imposes only two 
requirements: 

M 2 H - \M% = e Ac M| « Ml (31) 

and small Mx- The upper bound on the invariants involves both e\ and e c , and as we 
hinted above they partially cancel in the combination e\ c which enters into the require- 
ments: \e\ c \ < \e\\, \e c \. (One sum rule which can be formed under these boundary 
conditions, —2w^j + + |m^, is particularly sensitive to this cancellation 0.) The 
first requirement, Eq. (|3l|) , entails a definite tuning of parameters to a precision of 
€\ c , which does not apparently follow from any symmetry. Note that without this re- 
quirement, either color breaks when Mfj > \M^ q or a Coleman- Weinberg mechanism 
operates when < \M^ q . The requirement of small Mx may on the other hand be 
natural, since [see Eq. (fZ^f)] the value of M\ is smaller by an order of magnitude than 
the soft-breaking masses whose splitting generates the D-terms, and those masses 
may be expected to be comparable to M sq and Mh- In any case, we see that because 
of the focusing effect of the RG equations, the D-terms cannot be allowed to induce 
splittings bigger than those we already had through custodial SU(2)-breaking effects. 
Hence they do not eliminate the criticism that the electroweak symmetry is hard to 
break when the Yukawas are comparable. But there is still a significant advantage 
in using these D-terms, since they can now substitute for large values of m 1 / 2 and fi, 
and with light gauginos and fi it is much easier to obtain a large tan j3, according to 
Eqs. (|) and @T|). 



10 Radiative bottom decay (II) 

Before putting the various observations to use in examining specific scenarios and 
their merits, we point out another feature of the solutions to the RG equations which 



will further constrain the scenarios. As is evident from Eqs. (p~6| ~ p~8|) (or directly 
from the RG equations), the initial value Aq hardly affects the low-energy values of 
A ti b jT ; they are instead largely determined in magnitude and sign by the gaugino mass 
My 2 , which also fixes the low-energy gluino mass. (It is difficult, though perhaps not 
impossible for sufficiently small Ac, to construct models in which Ac ^> Mi/ 2 and 
yet the electroweak symmetry but neither color nor charge breaks spontaneously and 
correctly, so we shall disregard this possibility in these proceedings. The implications 
of tuning Aq to cancel the gaugino mass at low energies in the expression for A t will 
be considered elsewhere ||.) This observation, which was also emphasized by Carena 
et al. ||, directly relates the 5rrib corrections of Eq. (0) to the large b — ► graphs 
discussed above. (More precisely, the gluino- and higgsino-exchange diagrams for 
each process are directly related.) The sign of this correlation U is such that when 
Smi, < (i.e. the predicted R is lowered, and therefore so is the top mass) then the 
large b — ► 57 graphs interfere constructively with the usual 2-Higgs standard model 
amplitude, and vice- versa. On one hand, we see from Eq. (|3|) or from Fig. 1, that the 
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bounds on Srrib are more severe when Srrib > 0. On the other hand, as noted above, 
when Srrib < the interference is constructive and the bounds on the large b — > 57 
graphs are stricter. Thus these two bounds are much stronger when taken together, 
and translate into the following statement: either (a) the gauginos or fi or both are 
significantly lighter than the squarks, or (b) the superpartners are much heavier than 
the Z. 



11 Case studies 

With these remarks in mind, we first examine the popular |I] case of universal soft- 
breaking masses. This scenario has also been recently studied in some detail by Carena 
et al. ||. If all soft-breaking scalar masses are equal then the D-term contributions 
vanish (Mx = 0), and we are left with the three parameters /j,, M 1( / 2 and M = M sq = 
Mh, in addition to Bq which is adjusted at the end to obtain the correct tan/5 [see 
Eq. d§|)]. We have already mentioned that in this case we need /1 comparable to M 1 / 2 , 
and both at least as big as M , to break electroweak symmetry correctly. We have 
also seen that these three parameters must be tuned in order to obtain a positive m| 
and m\: 

m 2 z ~m 2 A ~ e c M 2 l/2 . (32) 

Next, to achieve a hierarchy of Higgs VEVs, Bq must be adjusted very precisely such 
that, at low energies, 

B 1 m\ e a 



(33) 

mi/2 tan p n mi/2 tan p 

Finally, since fi and the gauginos are not lighter than the squarks, Srrib is rather 
large, and so must be negative, as can be seen from Eq. (|3|) or from Fig. 1. Hence 
the b — ► 57 constraint is strong, necessitating large superpartner masses of at least 
0(TeV) and therefore a further tuning (by roughly another order of magnitude) of 
the three parameters to achieve correct electroweak breaking. Of course, such a 
highly-tuned scenario is also highly predictive: for example, the spectrum is highly 
constrained, and the top mass is predicted by the large Srrib corrections to be light. 

In search of a more natural scenario, we next relax the assumption of universal 
soft-breaking masses, which was perhaps arbitrary to begin with. We begin with 
the maximally symmetric scenario studied previously (see also scenario C above), in 
which n ~ M1/2 ~ A G ~ B G (~ m z ) <C m^, . . . ,m b (<g M sq ~ M H ) so that the VQ 
and 1Z symmetries are approximately obeyed. Since this hierarchy directly implies 
a small Srrib and therefore (from Fig. 1) a relatively large Ac the focusing effect of 
Eq. ( PB| ) takes its toll, and once again the initial parameters — this time M sq , M^, and 
M x — must be adjusted to at least O [max (eA, e c )]. Then, to truly get a maximally 
symmetric scenario, we choose to obtain a large tan /3 not by tuning to get a small B 
but rather by (equivalently) tuning to get a small m 2 z ~ 777^/ tan/3, which then allows 
small fi and 777i/ 2 relative to the typical low-energy soft-breaking masses and therefore 
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establishes approximate VQ and 7Z symmetries even at the electroweak scale: 

o 2 2 maX ( e A? e c) n/f 2 (oa\ 

- 2m v = m z ~ — M s (34) 

where Ms = Mh ~ M sq . After this adjustment, large tan/3 is automatic: 

B 1 ml 



-0(1). (35) 
mi/2 tanp/imi/2 

And no further adjustment is necessary to suppress b s'j since /x and Mi/ 2 are 
small. So this scenario requires much less adjusting than the universal case, although 
more than just the inevitable 1/ tan/5 tuning. (The actual tuning needed is in reality 
slightly more than indicated, due to the squark- and slepton-splittings induced by 
the D-terms and due to the effects of custodial symmetry violation, as mentioned 
above; qualitatively, though, the picture we described remains.) It is also predictive: 
the superspectrum is hierarchical with light charginos and neutralinos but heavy 
squarks, and since Xg is large so is the top mass. 

We can continuously retreat from this maximally symmetric case by increasing \i 
or Mi/2 ( or t ne related parameters), thereby losing VQ or 1Z, first at low energies 

1/2 

when /i or 7711/2 become comparable to e A Ms, and then at all energies when /j or 
mi/2 become comparable to Ms itself. Both may be interesting for model-building 
(for example, if \x is radiatively generated by A terms at the GUT scale, then VQ 
rather than TZ symmetry should be evident) and for comparison with experiments 
once the superspectrum is measured. In either case, the tuning is comparable to 
the maximally symmetric case, though less tuning is needed in Eq. ( p4|) and more 
in Eq. (|35|) . We will however defer their discussion to our more complete study [|J, 
and instead consider the case where both VQ and TZ symmetries are abandoned in 
favor of a smaller Xq- We will call this scenario, unimaginatively, the asymmetric 
case. It alleviates the focusing effect of Eq. (|26|) since now e A ~ 0(1), which in turn 
allows the initial conditions (specifically the D-terms) to split the Higgs doublets and 
the other multiplets by a large amount, so now Mx ~ Ms- Small Xq also entails 
larger (negative) bm^ corrections to correctly predict mi,/m r , and to this end we 
take \i ~ M1/2 ~ M sq ~ Mh = Ms- We then see, however, that what we gain by 
eliminating the focusing effect we lose by restoring the b — > 57 problem: since \x and 
mi/2 are 110 longer small, we are forced to raise the SUSY scale to ~ O(TeV), and 
therefore again to tune the initial parameters to make the Z light: 

- 2ml =m 2 z~ (t^Ms) ■ (36) 
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Since the D-term splitting of the Higgs is now large, m\ ~ M|, we find that Bq 
requires only the typical tuning 

_b_ = J_j^^J__ (37) 

777,1/2 tan (3 fi mi/2 tan/3 

This scenario is comparable in its naturalness (or lack thereof) to the symmetric case. 
The superspectrum is uniformly heavy rather than hierarchical, and the top is light 
since Xq is small. 
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12 Conclusions 



These last two scenarios are qualitatively the best one can hope for from an SO(10) 
model with Yukawa unification, for two reasons. First, obtaining a large tan/3 and 
thereby the top-bottom hierarchy is never natural in the MSSM, due to the LEP 
bounds on /i and m.1/2 which force either m 2 A to be much heavier than the Z or B 
to be much lighter than the gauginos. Second, the last two scenarios illustrate how, 
for large Xg, the inherent focusing property of the RG equations in the symmetric 
limit necessitates a further tuning of the initial parameters, while for small Xq a 
similar tuning is mandated by bounds on the rate of b — > 37. There are also in- 
termediate scenarios, with only one of the symmetries, but they are apparently no 
more natural. These various possibilities are considered in more detail elsewhere |§. 
The universal case is generically much more tuned than most of these scenarios, as 
we have shown. Hence departures from universality, although possibly dangerous for 
the flavor- changing neutral current interactions they can induce in some models, are 
strongly favored in achieving a large tan f3. 

What is the status of the predictions? Perhaps surprisingly, the top mass is not 
an independent prediction of Yukawa unification, but rather depends strongly (i.e. 
non-logarithmically) on a certain ratio of superpartner masses appearing in Smi,. 
Since the large- and small-Ac cases are equally fine-tuned, naturalness arguments 
do not single out any particular top mass within this SO(10) framework. Instead, 
information about the top mass can be combined with Fig. 1 and Eq. (|2|) to restrict the 
mass parameters which can be consistent with Yukawa unification, and to determine 
a favored superspectrum. It must be admitted that, with Yukawa unification, the 
attractive conventional picture of radiative electroweak symmetry breaking due to 
a large X t but small is largely lost: the symmetry is either broken radiatively 
using small custodial isospin-violating effects and extraordinarily fine-tuned initial 
conditions, or else the Higgs doublets are already split at the GUT scale. Furthermore, 
we have seen that all the large tan (3 scenarios are technically unnatural. On the other 
hand (and to some extent because of the necessary fine-tuning), they are certainly 
predictive and so will be tested in future accelerators: for example, if the charginos 
are light but the squarks are heavy, the top should also be heavy; if the SU(2)l- 
singlet bottom squark or the doublet sleptons were lighter than the other squarks 
and sleptons, this would be a sign that the D-term splittings were large [jlOf] ; and 
finally, tan/5 can itself eventually be measured to decisively confirm or dismiss the 
large tan (3 hypothesis. 

Even before any further experimental input, there are a few theoretical avenues 
worth pursuing which would make Yukawa unification much more attractive. First, 



as discussed also by J. Lykken in these proceedings ||TTJ, string models which lead to 



true grand-unified models as their low-energy effective Lagrangians can exhibit higher 
symmetries in certain sectors of the GUT than in other sectors. In particular, a string 
theory with an SO(10) gauge symmetry might break to an effective low-energy SU(5) 
GUT in a stringy way, leaving all three Yukawa couplings unified at the Planck scale 
a la SO(10) but splitting the 5 H from the 5 H soft masses, namely My from M|>. Note 
that unlike the D-term splittings we have considered before, and M|> could now 
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be split without necessarily decreasing some squark or slepton masses. This approach 
can thus provide more freedom in the choice of boundary conditions — although as 
we have shown, some features of the RG equations, in particular the focusing effect, 
are inherent in the equations themselves in certain limits and apply to any boundary 
conditions, so that freedom could either be the key to a more natural scenario (see the 
following remark) or could necessitate still more arbitrary fine-tunings at the GUT 
scale. Second, we have shown that in the symmetric case with large Xg, the ratio 
of the soft-breaking masses for the 10 H and the 16 3 needs to have a certain value 
if the electroweak symmetry is to break correctly. While such a value is arbitrary 
in the context of an SO(10) model and therefore apparently requires a fine-tuning, 
perhaps this value could be explained as a ratio of integral conformal weights in the 
context of a string theory into which the GUT is embedded. (Notice that this value 
is favored simply by the unification of Yukawas at some large scale, and does not 
depend on an SO(10) symmetry.) If such an explanation could be found, then the 
symmetric large-A^ case would now be strongly favored: it would only require the 
minimal 1/ tan/3 tuning (and would predict a heavy top!). In fact, if the squarks 
were to be experimentally determined to be heavy while the gauginos were light, then 
this case would be no more fine-tuned than the conventional small tan/3 scenario, but 
would have the advantage of explaining the top-bottom mass hierarchy (through the 
VQ and 1Z symmetries) — which, after all, was historically the motivation for studying 
Yukawa unification. 
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Figure Captions 

Fig. 1: The dependence of the low-energy values Xt,b,r an d of the ratio R = Xt>/X T 
on the initial condition Xf br = Xq, without any threshold corrections. The 
allowed range of -R exp t is shown shaded, and the minimal value of Xg allowed by 
this range in the absence of any corrections is indicated by the solid dot; lower 
values of Xq require finite, negative Srrib (see the text). The corresponding value 
of A t is marked by the shaded dot. The vertical scale on the right indicates 
the approximate tree-level top mass ~ 174A t GeV which would result from the 
values of A t on the left vertical scale; for example, the shaded dot predicts a 
heavy top, above 170 GeV or so. 
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Fig. 2: The leading (finite) 1-loop MSSM corrections to the bottom quark mass, 
namely Sm^. 

Fig. 3: Our predictions for the pole mass of the top quark, without superheavy 
corrections and using two qualitatively-different superpartner spectra, specifi- 
cally m higgsino ~ /j, = 100 GeV, m g i uino = 300 GeV, m wino = 100 GeV, m squark = 
^siepton = 1000 GeV and uia = 1000 GeV for which the Smf, corrections are 
small, and mhiggsino ~ /i = 250 GeV, m g i uino = 300 GeV, m wino = 100 GeV, 
m squark = ^-sicpton = 400 GeV and niA = 400 GeV, for which \Sm h /m b \ ~ 0.25. 
The upper or lower horizontal axes should be used for these two spectra, re- 
spectively. The "cloud" indicates the region where the theory becomes nonper- 
turbative at the GUT scale. Also shown are the estimated allowed mass ranges 
for the running parameter as extracted in our previous work 
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